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Abstract 

We show that certain free energy functionals that are not convex with respect to the usual 
convex structure on their domain of definition, are strictly convex in the sense of displacement 
convexity under a natural change of variables. We use this to show that in certain cases, the 
only critical points of these functionals are minimizers. This approach based on displacement 
convexity permits us to treat multicomponcnt systems as well as single component systems. 
The developments produce new examples of displacement convex functionals, and, in the multi- 
component setting, jointly displacement convex functionals. 

Mathematics Subject Classification Numbers: 49S05, 52A40, 82B26 

1 Introduction 

1.1 The variational problem 

We consider minimization problems for a type of functional that arises in the study of phase 
segregation in statistical mechanical systems. Let F{m) be a function on the real line that is 
continuous and strictly positive except at m = a and m = h with a <h. A. good example to bear 
in mind is the "double well potential" 

F{m) = ^(m^ - 1)2 , 

where of course a = — 1 and 6 = 1. 

Let Ca^b be the set of measurable functions m{x) from M to R such that (for some representative) 

lim m{x) = a and lim m{x) = h . 

x—*—oo x—>+oo 
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The numbers a and b represent the values of the order parameter m in two phases of a statistical 
mechanical system. For example, m = a might correspond to a vapor phase, and m = 6 to a liquid 
phase. 

A function m[x) in Ca^b denotes a possible one-dimensional transition profile across the boundary 
segregating the two different phases. The actual profile that one would expect to see would be one 
that minimizes the free energy cost of making such a transition. The free energy functional J- to 
be minimized on Ca^b will in some cases of interest have the form, c.f. |10j . 

T{m) = I F{m{x))dx -\ — / / {m{x) — m{y))'^J{x — y)dxdy , (1-1) 

where J{x) is a non-negative integrable function on M. 

The term F{m{x))dx is due to short range interactions and entropy effects and is normalized 
so that it vanishes in the pure phases, when m{x) = a or m{x) = b, while the term Jj^Jjg^(m{x) — 
m{y))'^J{x — y)dxdy is due to long range interactions. This long range term in the free energy 
suppresses sharp transitions, as does the gradient term in the familiar but purely phenomenological 
Van der Waals model |13j . For more discussion of the physical context of the problem, see [8]. 

Much useful information can be deduced from the specific form of the minimizing profiles. In 
particular, the surface tension at a two dimensional phase boundary in physical three dimensional 
space is the minimum value of J-{m) on Ca^b] see p] and Section [5] for more information. Hence we 
ask: 

• What is the minimum value of J-{m) as m ranges over Ca^bj (^nd are the minimizing profiles, if 
any, unique up to translation? 

Actually, the existence of minimizers is relatively simple to prove using the rearrangement 
inequalities to be discussed below. However, because of the translation invariance, they are never 
unique: Any translate of a minimizer is again a minimizer. It is less simple to show that this is the 
only degeneracy. 

1.2 Displacement convexity and uniqueness of fronts 

For a particular choice of F in the free energy functional specified in (|l.ip . the minimizing profile 
problem has been solved in a series of papers |10|.|llj by De Masi, Orlandi, Triolo and Presutti, 
building on previous unpublished work of Dal Passo and de Mottoni [9] Their solution involves the 
construction of a dynamics that is dissipative for the free energy functional, and then a careful 
analysis of limits along the time evolution for this dynamics. 

Another approach that we further develop here has been introduced by Alberti and Bellettini 
[2], [1]. They discovered an alternative convex structure which renders the variational problem for 
(jl.ip convex, and used this to study the existence problem in [2]. Later, Alberti [1] returned to the 
problem and proved a uniqueness result that affirmatively answers the question raised above for 
this one component model. 

Our goal here is to treat certain two component systems. Motivated by this problem, we 
were led to reconsider the single component problem from the point of view of McCann's notion 
of displacement convexity [12]. In fact, the minimization problem for (II. ip is challenging largely 
because the functional T is not convex on Ca^b in the usual way: For < A < 1, and mo and mi in 
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Ca^bi define m\ = (1 — A)mo + Xrrii and note that m\ S Ca^b- However, due to the non convexity of 
the potential function F, it is not true in general that T{m\) < (1 — A)^(mo) + AjF(mo). 

In [12], McCann, building on groundbreaking work of Brenier [6], introduced an alternative 
convex structure on the space of probability densities on M", and used this to prove existence 
and uniqueness results for minimizers of functionals that were not convex in the usual sense. We 
shall show here that the minimization problem for (II. ip . as well as for a two component model of 
this type, can be handled within this framework. In the process, we provide two new examples 
of strictly displacement convex functionals, the second of which is jointly displacement convex. It 
turns out that the alternative convex structure introduced in [2] is equivalent to the displacement 
convexity in this one dimensional setting, although the approach is quite different. We shall see 
that developing the alternative convex structure explicitly in terms of displacement convexity has 
advantages, especially for the two component system, when one seeks to prove a uniqueness result. 
Moreover, as we show in the final section, our results for the two component system may be 
applied to the single component system in higher dimensions, yielding a new uniqueness theorem 
for monotone solutions of the Euler-Lagrange equation. 

We now describe the alternate convex structure with respect to which T is convex. This second 
convex structure cannot be defined on all of Ca,b, but only on the subset A4a,6 consisting of right 
continuous monotone profiles. Nothing is lost in this restriction, as rearrangement inequalities 
show that minimizers of J- on Ca^b must actually be monotone, so that they have a right continuous 
version belonging Ma,b^ see [1] and Theorem lG.lI b elow . 

Any right-continuous profile m(x) in Ma^b can be written in the form 

m{x) = a+ {h - a) I d^(y) (1.2) 

where is a uniquely determined probability measure on R. This identification of M.a,b and the set 
of probability measures on M allows us to look at .F as a functional defined on probability measures. 

This is a useful perspective since there is an alternative convex structure on the set of probability 
measures on M (or more general domains) that was introduced by McCann, and which we describe 
below. A functional on probability measures is said to be displacement convex if it is convex with 
respect to this alternative structure. We shall show here that regarded as a functional on 
probability measures is, in fact, displacement convex. Using this, we shall show that any solution 
in Ma,b of the Euler-Lagrange equation for the variational problem concerning (jl.ip 

m{x) = — ( I J{x — y)m{y)dy — F' {m{x)) \ , (1-3) 

where 

J = / J{x)dx , (1.4) 
Jr 

is in fact a minimizer. Solutions to (jl.3p can easily be constructed by iteration and using these 
surface tensions may be readily computed. 

This solution to the variational problem has the advantage of applying also to free energy func- 
tionals in certain multicomponent systems, in which the determination of the minimizers has not 
been previously treated. Indeed, our motivation was to be able to rigorously determine the surface 
tension in such systems. However, we shall first present our simple solution of the minimization 
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problem for the single component free energy functional specified in (jl.ip . and then treat the 
multicomponent case. 



2 The alternative convex structure 
2.1 The reduction to monotone profiles 

First of all, notice that if we seek to minimize on Cq^^, we need only consider profiles m for which 
a < m{x) < b for all x. Indeed, for any m G C^^b, define fh by 

m(x) = min{6, max{a, m(x)}} . 

Then J^{fh) < J-{m) with equality only in case fh = m, since otherwise replacing m by m lowers 
both the potential and the interaction terms. 

We now recall a notion of rearrangement due to Alberti [T|. For any Borel measurable set A, 
let 1^1 denote its Lebesgue measure. The rearrangement is defined for Borel sets A C M such that 
|j4A(0, oo)| < oo, where A/S.B = A\B U B\A is the symmetric difference of A and B. For such a 
set A, define the rearranged set A* by 

A* = [a,oo) where a = |(0, oo)\^| - |^\(0, oo)| . 

Any function m in Ca^b that takes values in [a, h\ can be represented in "layer-cake" form: 



m[x) 



[ '^{m>z}{x)dz + a . 

J a 



For each z G (a, b), the set {m > z} certainly has the property that \{m > 2;}A(0, oo)| < oo. Hence 
one can define the rearrangement of m itself through m*{x) = J^^ {^{m>z})* {x)dz + a. (Applying 
the rearrangement to a monotone increasing function, one simply obtains the right continuous 
version.) 

Alberti shows that for any two such functions mi and m2, 

|m*(x) — m2(x)pdx < / |mi(j;) — m2(x)pdx . 

Jr 

In particular, with m being any function in Ca^b that takes values in [a, b], and h any real number, 
let mi{x) = m{x), and m2{x) = m{x + h). Then 

\m*{x) — m*{x + h)\'^dx < / \m{x) — m{x + h)\'^ dx , 

JR JR 

SO that 

\m*{x) — m* {x + h)\^dx J{h)dh < 



HI 

JR \JR 



\m{x) — m{x + /i)pdx ) J{h)dh . 

JR \JR 

This of course means that 

(m*(x) — m*(y))^ J(x — y)dxdy < / / {m{x) — m{y)Y J {x — y)dxdy . (2.1) 
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In fact, Alberti shows (see Theorem 2.11 in [T]) that there is equahty in (j2.ip if and only if m = m*. 

Of course, F{m*{x))dx = F{m{x))dx, and so we have T{m*) < T{m) with equality if and 
only if m = m* . Thus, we may restrict our search for minimizers to Ca,6, the subset of monotone 
increasing profiles in Ca^b- 

2.2 Displacement convexity of m i— j^F{m{x))(ix 

As noted in p.2p . if m is any profile in ^Aa,b■, then (m(x) — a)/{h — a) is the cumulative distribution 
function of a uniquely determined probabihty measure /u: 

m[x) — a 



o-a J(-oo,x] 

For each m in Ma,b^ define x(m) to be the inverse function: For m € (a, 6), 

x{m) = inf{ X : m{x) > m } . (2-2) 

Then of course, m{x) is the inverse function of x{m), so that for x in M, 

m(x) = inf{ m : x{m) > x } . (2-3) 

Let dx(m) denote the Lebesgue-Stieltjes measure on [a, h] induced by the monotone function x{m). 
(In the terminology introduced below, dx(m) is the push- forward of Lebesgue measure on M under 
m.) Then one can rewrite 

.b 

F{m{x))dx = / F(m)dx(m) . 



Let mo and mi be any two elements of A^a.fe, and let xq and xi denote their respective inverse 
functions. Then for any A G (0, 1), define x\{m) by 

x\{m) = (1 — A)xo(m,) + Axi(m) . (2-4) 

Note that xx is also the inverse function of an element of Ma,b-, which we shall call mx. That is, 

mx{x) = inf{ m : (1 — A)xo(m) + \xi{m) > x } . (2-5) 

Note that dxx{m)^ the the Lebesgue-Stieltjes measure on [0, 1] induced by the monotone function 
xx{m)^ satisfies dxxirn) = (1 — A)dxo -|- Adxi. Then, 

F{mx{x))(lx = I F{m)(lxx{'m) 



rb rb 

(1 — A) / F{m)dxQ{m) + A / F(m)dxi(m) 



(1-A) / F{mQ{x))dx + \ / F{mi{x))dx . 

(2.6) 



This tells us that along the interpolation mx between mo and mi provided by (j2.5p . the function 
A 1-^ ji'^F{mx{x))dx is affine, and in particular, is convex. This is not the case for the standard 
interpolation given by 

mx{x) = (1 — A)mo(x) + \mi{x) , (2-7) 



CCELM July 2008 



6 



since A i-^' J^F{fhxix))dx is not, in general, convex. That is, taking convex combinations in 
terms of the inverse function x{m), as in ()2.5p . instead of m{x) itself, as in (j2.7p . has "cured" the 
non-convexity of the functional m i— > J^F{m{x))dx. 
Of course, this will only be useful if the functional 

»«//(™W-™te))V(.-.)d.d., (2.8) 

Jr Jr 

which was convex in the usual way, is still convex with the new convex structure. This is not at all 
obvious, but the main result of the next section asserts that this is the case. 

The approach of Alberti and Bellettini [2], which we discovered only after our work was complete, 
was to rewrite the interaction directly in terms of Xm, and to show that it is convex. 

However, it turns out that the convex structure in (12. 4|) is something that is by now well- 
known; it is the displacement convexity structure introduced by McCann. Making this connection 
will facilitate showing the strict convexity of m i-^ /k /m("^('^) ~ in(y))'^ J{x — y)dxdy under this 
convex structure. This point was left open in [2], who explicitly asked whether one could extend 
the ideas to give a direct proof of uniqueness. Although Alberti [T| did later return to address the 
issue, we shall see here that the strict convexity is quite clear from the perspective of displacement 
convexity. 

Displacement convexity is usually introduced as a convex structure in a set of probability mea- 
sures. Given a probability measure /Uq on M, and a measurable map T : R — > M, we define the push 
forward of fio under T, T^fiQ, by 

4>{T{x))dMx) = [ 4>{y)d{T#f^o){y) , (2.9) 



for all bounded, continuous functions (p. 

Given two probability measures /xq and jii on M, there is a unique monotone map T such that 
T^fiQ = To see what it must be, fix any a G R, let (j)a be the step function cpaix) = l{-oo,a] (2^)- 

Then, by definition, we must have / 4>a{T{x))dfio{x) = / (j)a{y)dfj,i{y), and hence 

Jm. Jm. 

/ d/xo = / dA^i . (2.10) 

J —00 J —00 

Let mo and mi be the cumulative distribution functions of /io and ^1, respectively. Then (j2.1U|) 
entails that mo(T~^(a)) = mi (a) for all a, or, what is the same thing 

mo(a) = mi(r(a)) (2.11) 

for all o. As long as mi is free of "flat spots" , so that the inverse function does the expected thing, 
this leads to 

r(a) =xi(mo(a)) . (2.12) 

As long as ^iq and /ii have strictly positive densities, (|2.12|) does indeed define a monotone map T, 
and then it is very easy to see that with T defined by (|2.12|) . T^fiQ = /ii, and in fact, this is true 
without further technical hypotheses; see [13] for more information. 
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We now interpolate the map T, and hence the corresponding probabihty measures and /xi 
and the corresponding cumulative distribution functions tuq and m-i as well. For all A G [0,1], 
define Tx by 

Txix) = (1 - X)x + XT{x) . (2.13) 

If we define x\{m) by 

x\{m) = Tx{xo{m)) , 

then clearly xx is given by (|2.4|) . 

The displacement convex structure on probability measures on M is given by fix = Txi^fJ-o, and 
so it is nothing other than the convex structure (|2.4p . expressed in terms of probability measures 
instead of cumulative distribution functions. When /xi and /i2 have strictly positive densities, so 
that T is given by (j2.1ip . we denote the density of fix by px, and write 

px = Tx#Pi ■ (2.14) 
We summarize the main result of this section in a theorem: 

2.1 THEOREM. Let A i-^ mx be the displacement interpolation between rriQ and mi in Aia,b- 
Then /or < A < 1, 

-F(m;^(x))dx = (1 — A) / F{mQ{x))Ax + \ / F{mi{x))dx . 

Jr Jr 

3 Displacement convexity of the interaction energy 

Let M denote the class of cumulative distribution functions on M. Making the obvious change of 
variables, we will assume without loss of generality that a = and 6 = 1 and we will set Mo,i = 

Given any m G A1, let /u be the corresponding probability measure, so that m(x) = dfi{y). 
The first step in the investigation of the interaction energy is to rewrite it as a functional of fi 
instead of m. This is done in the following lemma: 

3.1 LEMMA. Assume that \s\J{s)ds < oo. Define W in terms of J by setting 

/•oo 

W{u)= (s-n)(J(s) + J(-s))ds . (3.1) 

J u 

for u >0 and W{u) = W{—u) for u < 0. Then 

{m{x) — m{y))^ J {x — y)dxdy = / / W {z — w)dfj,{z)dfi{w) . 
I Jr Jr 

W is a symmetric function, and is convex on (0, oo) and on (— oo,0), though not on all o/M. 
Proof: Since for x < y, m{x) — m{y) = JJd/i(z), we have from the Fubini Theorem that 

(m(x) — m{y))'^J{x — y)dxdy = 



^[x,y]iz)l[^,y]{w)dp{z)dp{w) 



J{x — y)dxdy 



'^[x,y]iz)l[^,y]{w)J{x - y)dxdy 



dfi{z)dfj,{w) . 

(3.2) 
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Thus if we define V{z — w) by V{z — w) = / l^x,y]{^)^lx,y]{'w)J{x ~ y)dxdy we have 

{■m{x) — m{y))'^ J{x — y)dxdy = / / V{z — w)dfi{z)dfj,{w) . 
I Jr Jr 

We next show that V = W. To do this, write 

) J{x) for X > , 
J+[x) = < and J_(x) = J(x) — J+ix) . 

[O for X < 

We first consider l[x,y]{z)l[^^yj {w)J^{x — y)dxdy. Make the change of variables s = y — x and 
t = {x + y)/2. Then dxdy = dsdt, and l[^^y]{z)\^^y^{w) = l[t-s/2,t+s/2]{z)l[t-s/2,t+s/2]{w). This 
quantity is zero unless |z — tf] < s and \2t — {x + w)\ < s — \z — w\, in which case it is one. Therefore 

/■oo / i-{z+w)/2+{s-\z-w\)/2 \ 
/ dt \ J+{s)ds (3.3) 

,z-w\ \J {z+w)/2-{s-\z-w\)/2 I 



oo 



{s — \z — w\)Jj^{s)ds 



Doing the same calculation for the part involving J_, we obtain that V = W where W is 
given by ([ST]). Also, for n > 0, W'{u) = - J!^{J{s) + J(-s))ds, and so W"{u) = J{u) + J{-u), 
which is non-negative. Thus, W is convex on (0, oo), and on (— oo,0) by symmetry. However, 
it is not convex on the whole real line. Notice that W{Q) = s{J{s) + J{—s))ds > 0, while 
lim„^±oo W{u) = 0. □ 

We now prove the main result of this section: 



3.2 THEOREM. Let A i— > m\ be the displacement interpolation between mo and mi in A4, as 
defined in h2. 5|) . Then for < A < 1, 

{mx{x) -mx{y)fJ{x -y)dxdy < (1 - A) / l {mQ{x) - mo{y)f J {x - y)dxdy 

JR JR 

+ ^ {mi{x) — mi{y))'^J{x — y)dxdy . 



(3.4) 

// J is strictly positive on some interval, and mo has a strictly positive derivative almost everywhere, 
there is equality if and only if mi is a translate of mo. 

Proof: If W were convex on all of M, the displacement convexity of the interaction energy would 
be a classical result of McCann [1^. However, in one dimension, the partial convexity of W that 
was established in Lemma [3TT] suffices, as observed by Blower |5|. This is because the map Tx is 
monotone for all A. Therefore, if z > w, Tx{z) > Tx{w) for all A. Hence, as we vary A, Tx{z) — Tx{w) 
stays in a domain of convexity of W . 

Hence from (|2.4p . if djjx = Txi^dfiQ is the displacement interpolation between d^o and d/ii. 



// 

JR JR 



Wiz-w)dfixiz)df,xiw)= WiTxiz)-Txiw))d^loiz)d^,oiw) . (3.5) 
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Define the map S{x) by S{x) = T{x) — x. Then, we can rewrite (j3.5p as 

W{z - w)dfix{z)dfix{w) =11 W{[z -w\ + \[S{z) - S{w)])d^,o{zWo{w) . (3.6) 



By the remarks made above, the right hand side is clearly a convex function of A. In fact, under 
mild assumptions on /^o or J, it is strictly convex unless T is simply a translation. 

To see this formally, let J be symmetric for simplicity of notation, and differentiate the right 
hand side of (j3.6p twice in A, finding 

2J {[z -w]+ \[S{z) - S{w)]) [S{z) - S{w)fd^iQ{z)d^lQ{w) . 



If this vanishes for all A, then / / 2J{z — w)[S{z) — S{w)]'^dfio{z)dfiQ{w) =0. If J is strictly 

Jm. Jm. 

positive and if /^o has a strictly positive density, then this is possible if and only if S is constant, 
and that of course means that T is a translation. 

To make this argument rigorous, and to relax the hypotheses, let /(A) denote the right hand 
side of (|3.4p minus the left hand side. Then, with g{z,w,X) defined by 

g{z, w, A) = [XW{z -w) + {l- X)W{{z - w) + {S{z) - S{w)))] 

-Wiiz - w) + XiSiz) - S{w))) , 

we have /(A) = J^J^g{z,w, X)dfio{z)dfio{'w). Since the integrand is non negative, we have for any 
measurable subsets A and B of M, 

/(A)> / [ g{z,w,X)dMz)disoiw) . (3.7) 
J A Jb 

Suppose that J is strictly positive on the open interval I = {i/q — 5/2,yQ + 5/2). Then / is an 
interval of strict convexity of W, so that whenever z/; — z e /, A i— > g{z,w, A) > on (0, 1) unless 
S{z) = S{w). However, if S is not constant almost everywhere, we can find an arbitrarily small 
interval about some zq on which it has strictly positive oscillation. In particular, we can find a 
zq and an e > so that /^^|'J*'//''2^(5'(z) — c)'^dz > e for all c. Let A = {zq — 6/2, zq + 6/2), and let 
B = {yo + xo — 6/2,yo + xo + 6/2). Then for all z in A and w in B, z — w belongs to /. Moreover, for 
every w in B, f^{S{z) — S{w))'^dw > 0, so \S{z) — S{w)\ > on a subset of A of positive Lebesgue 
measure. Since no has a strictly positive density, this ensures that the right hand side of (|3.7p is 
strictly positive. □ 

It is clear that the conditions on J and mo that are invoked to ensure strict convexity can be 
relaxed, though they are already quite general. 

We close this section a remark. If the profile m is continuously differentiable with m'{x) = p{x), 
and J{x)dx = 1, then 

lim / / -K — J I — r — dxdy = p [x)dx . 



h^oJ^J^ /l2 h \ h ^ 

It is already well known that the functional p p'^(x)dx is displacement convex, so the fact that 
Theorem 13.21 gives another proof of this is not of great interest. However, the connection between 
the two functionals at least gives one suggestion as to why the interaction functional might be 
expected to be displacement convex. 
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4 For the functional JT, critical points are minimizers 

4.1 THEOREM, //mo is any critical point of J- in M, and m is any other profile in A4, then 
T{m) > T{mQ) and there is equality if and only if m is a translate of ttlq. 

Proof: Let m\ be the displacement interpolation between tuq and m. Then A ^ J^{m\) is convex, 
and the derivative is zero at A = 0. Hence tuq is a minimizer of J-^ so that J-{m\) > J^{mo), and 
there is equality if and only if A i-^ J-{m\) is constant. But in this case, the strict displacement 
convexity of J- ensures that m is a translate of mo. □ 

5 Fronts in a binary fluid model 

We now turn to the study of the analogous problem for a binary fluid model. The binary fluid 
model has been investigated in [7] and [8], and we refer to those papers for details. Although the 
arguments apply to that setting in full generality, we discuss here only a special case where the non 
local interaction is only between particles of different species and the local term is purely entropic, 
for the sake of brevity. For further information and a numerical investigation of the minimizing 
fronts, see pj. 

In what follows, m(x) and n{x) represent the particle number densities of two different species 
of particles contained in some bounded domain J7 in M". Consider the functional J- defined by 



Here, J is a non negative, decreasing and compactly supported function on with range R. 
Notice that we must impose more conditions on J in the case of two species that we did in the 
single component model. The reasons for this will be made clear in Section 6. 

The problem considered in [7] is to minimize J-{m, n) subject to the constraint that 



have certain prescribed values. As shown in [7], this system undergoes a segregating phase transition 
when beta is large enough for the interaction term to overcome the entropy terms in T. These would 
prefer to have m and n to be uniform and this will indeed be the minimizing state for small /3, i.e. 
high temperature (3~^. However, for large values of /3, the advantages of segregation can dominate, 
and the fluid separates into two phases, one rich in particles of species 1, and the other rich in 
particles of species 2. Our concern here is with the proflles of the densities at the interface between 
the two phases. 

The nature of the two phases in the bulk is determined by considering the zero range model, in 
which the length scale R of the interaction J is negligible compared to the size of O. Formally this 
corresponds to setting J(x — y) = J6{x — y). It is also convenient to drop the constraint (|5.2p and 
to consider the function 






and 





//3,Ai,A2 ("^) "^j = m In m + n In n + (3Jmn — Aim — 



(5.3) 
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as a local free energy density. Here, as in the one component case, J = f^n J(|x|)dx , and Ai and A2 
are Lagrange multipliers that ensure the constraint (j5.2p on the total particle numbers. One may 
also think of Ai and A2 as specified chemical potentials and then determine m and n as functions 
of Ai and A2. 

In |7j it is proved that, if Ai 7^ A2, there is an unique couple fh,n minimizing fp^x^^x^ifn^n). 
However, if Ai = A2, there is a f3c such that, if P < Pc the minimizer is still unique, while, \i (3 > f3c 
there are densities < such that the couples {p'^,p^) and {p^,p^) are both minimizers of 
//3,Ai,A2 ("^) ''^)- We focus on the last case. Thus, in what follows Ai = A2 = A. Then fp^XiM is the 
local Gibbs free energy density. 

Analysis of the zero range model suffices to determine the quantity of the fiuid that is present 
in each phase, but not the surface tension across the boundary. We now turn to the variational 
problem that determines the density profiles across the interface, and the surface tension. We will 
assume that the geometry of is such that the interface is perpendicular to the first coordinate 
axis; e.g., we take to be a very long cylinder along the xi-axis with periodic boundary conditions 
along the other coordinate axes. 

First, we need to introduce the one dimensional version of J. Choose coordinate (s,t) on M" 
with s G M and t G M"'"^ and define J on R by 

J[S)= I J(^s2 + |t|2)dt , 

and then J = J{s)ds. Let g/^^x = infm,n>o //3,a,a("^7 By what has been noted above, 

5/3, A = ft3,x,xip',P^) = ff3,x,xip^,P~) ■ 
The functional G defined by 



G{m, n) 



m{x)\nm{x) + n{x)\n.n{x) + (3 I J{x — y)m{x)n{y)dy — gj3^x 



dx (5.4) 



is the excess free energy of a front. We look for the minimizers of this functional for /? > The 
minimum value gives the surface tension across the planar phase boundary. Note that we have let 
— > M and that G is the free energy per unit {d — 1)- dimensional area. 

Our goal in the next sections is to prove a strict displacement convexity property of this excess 
free energy functional, and to show, as a consequence, the uniqueness of the minimizing fronts up 
to translation. As in the one component case, a rearrangement inequality will enable us to restrict 
our attention to monotone profiles. Let ^Ap~ p+ be the subset of Cp- p+ consisting of monotone 
increasing profiles, let Aip+ p- be the subset of Cp+ p- consisting of monotone decreasing profiles 

Our main goal mathematically in what follows is to show that the functional 



(m, n) 



J{x — y)m{x)n{y)dy — Jp^ p 



dx 



is displacement convex on Aip- ^p+ x Mp+ p- , where now we have both an increasing and a decreasing 
density profile. 

We shall prove the displacement convexity results in the next section. This time, we shall require 
certain moment conditions to obtain the displacement convexity. Hence, before we can apply these 
results, we need to show a priori that all minimizers have good localization properties. We do this 
by an analysis of the Euler Lagrange equation. 
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5.1 Convexity of the interaction energy for Q 

Define the functional I on Aia b ^-^crfby 



X(m, n) 



dx 



J{x — y)mi{x)m2{y)dy — jrh{x)n{x) 



R 



(5.5) 



We assume J to be non negative, even and compactly supported on M with range R. We define J 
to be the total mass of J, and we define 



m(x) 



b for X > , 
a for X < 



and 



n X 



d for X > , 
c for X < . 



Note that in the special case a = d = p and b = c = p^ , 



Z{m, n) 



J{x — y)m{x)n{y)dy — Jp^ p 



dx 



Although this special case is all that is needed for our applications here, we treat the general case 
because the small extra efi^ort yields a broad new class of jointly displacement convex functionals. 

The first step in our analysis is to rewrite I as a functional on probability densities. Let the 
probability densities pi and p2 be defined by 



m(x) = a + (6 — a) / pi{t)dt; n(x) 



c+{d-c) I P2{t)dt 

J — oo 



(5.6) 



We the rewrite the functional in terms of p\ and p2, and integrate by parts. Formally, one moves 
an antiderivative from each of pi and p2 over to J. Since J is positive, integrating it twice produces 
a convex function W ^ different from the one constructed in the one-component case. This is indeed 
what happens, but one must be careful about the boundary terms. The boundary terms do not 
vanish, but as we shall see, they depend on the densities in a very nice way, and altogether, one 
obtains the desired displacement convexity. 
To carry out this analysis, define 

\ f { I J{s)d.^ dt for X > , 
W{x) = <^ io VVo y (5.7) 

W{-x) for X < . 

Then W"{x) = J(x), H^(0) = 0, and W is an even convex function. Furthermore, 

lim W'ix) = — , Wix) = a + —\x\ for [x| > i? . 
x^oo 2 2 



(5.8) 



5.1 LEMMA. Let m € Ma,b o-nd n E Mc,d- Let pi and p2 be the corresponding probability 
densities defined in \5. &^) . Then, provided pi and p2 have finite first moments, and with W and a 
defined as above. 



Z{mi,m2) = {a — b){d — c 

+ [2{b - a){d - c) + be + ad]a 
J 



W{x - y)pi{x)p2{y)dxdy 
'. + ad]a 

X [(6 + a){d — c)p2{x) + (6 — a)(c + d)pi{x)\ dx 



(5.9) 
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Note that {a — b){d — c) > for b > a and c > d, which is the case when a = d = p~ and 
c = b = p^. Thus, (a — b){d — c)W{z) is a convex function of z on all of R. It follows in the 
usual way that the first term on the right is displacement convex. Since W is strictly convex on 
the support of J, it follows as in the proof of Theorem 13.21 that this part of the functional (j5.9p is 
in fact strictly convex apart from translation. The second term on the right of (j5.9p is a constant. 
The third term is a linear combination of the first moments of pi and p2- Since these first moments 
are displacement affine, we see that altogether, I{m, n) is strictly displacement convex, apart from 
translation. 

The fact that Lemma |5 . 1 1 requires a conditions on first moments, while Theorem 13.21 does not, 
means that it will be a little more work to apply Lemma I5.lt We shall need an a priori estimate 
guaranteeing that for any critical point (m, n) of Q, the corresponding densities have finite first 
moments. We shall return to this after first proving the theorem. 

Proof: We start by considering the integral in x first, on a bounded interval [—L,L]. Since 

J{x — y) = — W{x — y) we have that 
oxoy 



Moreover, 



dxdy 



W{x — y)m{x)dx 



-^W{L-y)m{L) + ^W{-L 



-L dy 
y)m{-L) 



y){b - a)pi{x)(ix 



(5.10) 



L rL 



J{x — y)m{x)n(y)dy 



LJ-L 
L 

L 



L rL Q 

—W{x — y){b — a)pi{x)n{y)dxdy + 



LJ-L 



dy 



■— T^(L - y)m{L) + q^W{-L - y)m{-L) 



n{y)dy 



(5.11) 



Now we integrate by parts once more, this time in y: 



^ d 

-^W{x - y)n{y)dy 



L dy j-L 
+W{x - L)n{L) - W{x + L)n{-L) . 



W{x - y){d - c)p2{y)dy 



(5.12) 
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Summarizing, 



L rL 



J{x- y)m{x)n{y)dy 



L J-L 



= -{b-a){d-c) / W{x - y)pi{x)p2iy)dxdy 

J-L J-L 

+{b-a) J [W{x-L)n{L)-W{x + L)n{-L)]pi{x)d 



— VF(L - y)m(L) + " y)m{-L) n(y)dy 



(5.13) 



Let us examine the boundary terms 

r-L 



-L [ dy 



W{L - y)m{L) + ^W{-L - y)m{-L) 



dy 



n{y)dy , 



We have 



B2 := (b-a) J [W{x - L)n{L) - W{x + L)n{-L)] pi{x)dx 
Bi = m{L) W{L -y){d- c)p2{y)dy - m{-L) W{-L -y){d- c)p2{y)dy 

J-L J-L 

+m{L) [-W{L - y)n(y)]+i + m{-L) [Wi-L y)n{y)]ti • 
= (d - c) [m{L)W{L -y)- m{-L)W{-L - y)] p2{y)dy 

+m(L) [-W{0)n{L) + W{2L)n{-L)] + m(-L) [W(2L)n{L) - W{0)n{-L)] 
For 2L > R, where R is the range of the interaction J, the last two terms give 

{bc + ad){JL + a) + 0{l) 

To compute the other term, we consider, for a function / rapidly decaying, Jj^^ f{x)W{x + L)dx 
and J^j^ f{x)W{x — L)dx. We have 

rL r-2L r-R j 

I f{x)W{x + L)dx= I f{z-L)W{z)dz= / f {z - L)W {z)dz + / f{z-L){-z + a)dz 
J-L Jo Jo Jr 2 

The first term vanishes in the limit L — > 00 because of the decay of / and of the boundedness of 
W{z) for z G [0, i?]. The second term becomes, if j'^ |x|/(x)da; < cxd, 

r f{^)d{x + L) + a)dx='^ [ xf{x)dx + {a + ^L) [ /(x)dx + 0(l) 
J R-L ^ ^ ^ Jr 

In conclusion, 

/ j{x)W{x±L)dx = ±.'^- [ xf{x)dx + ia + L'-^) [ f{x)dx + 0{l) 
J-L ^ Jr 2 7r 
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Now we apply this result to B2, where the decaying function is pi, to get 
B2 = {b-a 
Now we apply to Bi: 



(c + d)— j X pi{x)(ix + a{d — c) f pi{x)dx + —L{d — c) j pi{x)dx 
2 ./TO ./to 2 



+ 0(1) 



Bi = {d- c) 



+ / xp2{x)dx+ (5.14) 

+ {bc + ad){JL + a) + 0{l) 



a{b-a) I p2{x)dx + ^L(b - a) j p2{x)dx 



Finally, 



Bi + B2 — J / m{x)n{x)dx = [2(6 — a){d — c) + be + ad]a 
Jr 

^ib + a){d-c) [ yp2{y)dy-i-{b-a){c + d) [ xpi{x)dx + 0(1) 

^ ./» ^ ./IB 



□ 



Lemma 15.11 is the key ingredient to prove the analog of Theorem 13.21 for the two-component 
model introduced in the beginning of this section. We now return to this model, and shall apply 
the lemma with a = d = p~ and b = c = p'^. Let {wi,W2) and (^1,^2) be in A4p- p+ x M^+ p- , with 
corresponding probability densities (r/i,r/2) and (C11C2), and let Ti, T2 be the monotone maps such 
that (i = Ti^rji, i = 1,2. Moreover, let A 1— > {rnxjUx) be the displacement interpolations between 
{wi, W2) and (^1,^2) and Tl{x) = Ax + (1 — \)Ti{x). 

5.2 THEOREM. Suppose that the probability densities iji and Q, i = 1,2 have finite first moments. 
Then for < X < 1, 



g{mx,nx) < (1 - A)^(u;i, zi;2) + A^(ui, ■U2) • 



(5.15) 



// J is strictly positive on some interval, and {wi,W2) have strictly positive derivatives almost 
everywhere, there is equality if and only if {vi,V2) is a translate of [wi,W2)- 



Proof. Lemma |5. II is applicable by the assumption that the probability densities have finite first 
moments. We set Si{x) = Ti{x) — x, so that 

l{mx,nx) = {p+ - p-f I [ W[ix-y) + X{S,{x)-S2{y))]drjix)dr,{y)+A{mx,nx) , 

Jr Jr 

with A affine. The function W is convex on all M, thus the interaction part of the Q is strictly 
displacement convex. Then remaining term is simply a linear combination of functions of m and n 
to which we can apply Theorem l2.11 The strict convexity up to translations follows as in the proof 
of Theorem O □ 
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Remark: In the two component case we need to use two monotone maps instead of one as in 
Theorem 13.21 Therefore it is crucial that W is convex on all of M and not just on (0, +oo) and 
(—00,0) as in the one component case. 

We close this section with a corollary showing that one could also use Lemma 15.11 to prove 
displacement convexity of the interaction energy in the one component model. In fact, in this 
application, the first moment condition drops out. 



5.3 COROLLARY. Let J satisfy the conditions below 115. 5\) . and and W defined as in the ^5.7\ ) 

Let m he a function that increases monotonically from —mp to m/^. Let p denote m' , i 
of m. Consider the functional <I>(m) given by <I>(m) = J{x — y) m{x)m{y) — 

<I>(m) = — 4m^ / / W{x — y) p{x) p{y)dxdy — Qarnjj . 



Proof. The functional ^{m) is equal to —Q{mi,m2) by putting mi{x) = m{x) and m2{x) = 
—m{x). This shows that — $ is strictly displacement convex, up to translation. □ 



6 Properties of the minimizers of Q. 

We restrict our attention to the case a = d, b = c. We need two results on the minimizers for Q, 

the first of which allows us to restrict our attention to monotone profiles when seeking to minimize 
Q. The second guarantees the existence of moments for the two densities corresponding to any 
minimizing pair (m,n). These theorems are: 

6.1 THEOREM. Suppose that J{x) is even non negative and decreasing. Then any minimizer 
(mi, 7712) of Q {mi, 1112) in Cp- p+ x Cp+ p- is monotone in the sense that mi is increasing and m2 
is decreasing. 

This theorem makes it easy to establish the existence of minimizers for G. The minimizers 
satisfy an Euler-Lagrange equation from which we can deduce a priori estimated needed to apply 
Lemma 15.11 

6.2 THEOREM. Suppose that J{x) is even non negative and decreasing on M+. Any minimizer 
w = {wi,W2) of Q in Cp- p+ x Cp+ p- satisfies p~ < Wi{x) < p+ for any x € M. It has derivative 
almost everywhere which is strictly positive and with ||u'^||li(k) is bounded. Furthermore, it satisfies 
the Euler-Lagrange equations 

lnm{x) + f3{J * n){x) = p, lnn{x) + (3{J * m){x) = fi, (6.1) 

where p = fii — 1 and * denotes convolution. Its derivative w satisfies almost everywhere the 
equations 

<^ + PiJ*w'2){x)=0, !^ + /?(J*^)(x)=0 (6.2) 

Finally, it converges to its asymptotic values exponentially fast, in the sense that there is a > 
such that {wi{x) — p^)e"l^l ^0 as x — > ^00 and {w2{x) — p±)e"^^^ — > as x — > =Foo. 
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The proof of Theorem l6.1l is adapted from a related result in for functions on the d dimensional 
torus. One could instead adapt the proof of Alberti's rearrangement inequality in [1] and remove 
the requirement that J be decreasing. But the present approach has the advantage of working also 
on the torus, and not only the line. The proof of the final part of Theorem 16.21 which is important 
for our application here since it provides the existence of moments, is adapted from the proof of a 
similar result for the one component system in pTj. In the rest of this section, we present these 
proofs. 

Proof of Theorem I6.lt To show this, we use a rearrangement inequality similar to those intro- 
duced in [7j for the analogous problem in the d-dimensional torus. For any xq € M, let T^^ denote 
the reflection about xq: 

Txq{x) = 2xq-x . 

Then define P, as the set of functions on M having finite limits at ±00 and the operators R^^ on 
Pby 

.^ /max{5r(x),5r(r^J} if x > xq , 
R^^g{x) = < (6.3) 
[mm{g{x),g{T^g)} if x < xq. 

c- r / N /max{/i(x),/i(T,J} if x < xq , 
[min{/i(x), h{Tx^)} if x > xq. 

lim g{x) if X < 
Let us also define ^(x) = aj^j fi similarly. 

lim g{x) if X > 

For any fixed xq and g,h & V, let g* denote R^^^g and h^, = R^^h. We now wish to show that 



g{x)J{x - y)h{y)dy - Jg{x)h{x) 



dx > 



g*{x)J{x - y)K{y)dy - Jg{x)h{x) 



dx 



with equality if and only \i g = g* and h = h*. 

To do this, let denote the half line {x [ x > xq}, and EI_ denote the half line {x | x < xq} 
and observe that 



// 



g{x)J (x - y)h{y)dxdy = 
g{x)J{x - y)h{y)dxdy + I I g{x)J{x - y)h{y)dxdy + 



g{x)J{x-y)h{y)dxdy+ I I g{x) J {x - y)h{y)dxdy = 
g{x)J{x -y)h(y)dxdy + / g{T^^,x) J {x - y)h{T^ay)dxdy + 
g{TxoX)J{T^^x -y)h{y)dxdy + I I g{x) J [x - T^^^y)h{T^^y)dxdy 



(6.5) 



The desired inequality is then a consequence of the following inequality for pairs of real numbers: 
Let oi and 02 and 61 and 62 be any four positive real numbers. Rearrange ai and 02 to decrease, 
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and 61 and &2 to increase; i.e., let a\ = max{ai,a2}, 02 = min{ai,02}, b\ = min{6i,62} and 
62 = max{6i,62}- Then 

+ a-^b^ - aibi - 0262 = A < 0, (6.6) 

462 + ^2^1 - «i&2 - 02^1 = -A > 0, (6.7) 

and there is equahty if and only if ai = and 61 = b* or a\ = 02 and 6^ = 62- 
We now apply the above inequalities with 

ai = g{x) 02 = g{TxQx) bi = h{y) and 62 = K^xoV) ■ (6.8) 

Then 

a\ = Rt^g{x) a*^ = Rt,g{T,,x) b\ = R^^h{y) and b''^ = R~^h{T,,y) . (6.9) 
Since J(r^(,x -y) = J{x - T^^y) < J{x - y), we get 

g{x)J{x - y)h{y) +g{T^ox)J{x - y)h{T^Qy) + 
9{Txox)J{T^^^x - y)h{y) +g{x)J{T^„x - y)h{T^^y) - 
Rto9{x)J{x - y)R^^^h{T^^y) + Rt^g{T^^x)J{x - y)R~^h{T^^y) - 
Rto9{Txox)J{T^^,x - y)R^^^h{y) + RX,g{x) J {T^^x - y)R^^^h{Tn,^,y) = 

-^[j[x-y)-J{x-T^,y)\> 

(6.10) 

for almost every x and y in IHI^-, with equality if and only if 

g{Txox) < g{x) and KT^^y) > h{y) (6.11) 

or 

g{Tx,x)>g{x) and h{T^,y) < h{y) (6.12) 

for almost every x and y in M^. Now unless g is constant, we can find x and xq so that either 
g{TxQx) < g{x) or g{TxQx) > g{x). Suppose it is the first case. Then (|6.1ip holds, and for almost 
every y, we must have h{TxQy) > g{y)- Making a similar argument for we see that one of (j6.1ip 
or (I6.12P must hold for almost every x and y. The only way that this can happen is if g and h are 
monotone. Now, by integrating (j6.10p on we conclude the proof. □ 

Proof of Theorem 16. 2t Everything but the exponential decay is standard, and details of the 
proofs of similar results can be found in [7J. To prove the exponential decay, we once again take 
advantage of the finite range R of J. 

Define a transformation $ : ^ by $(m,n) = (^Q^J■-^^Jn ^ ^fi-i3Jmy Then {p^,p~) and 
(p~ , p'^) are two stable fixed points of 'I'; the Jacobian of <I>, D^, is a strict contraction at either of 
them. Thus, there is a 5 > and an e > so that if 

\m — p'^l + \n — p~\ < 6 =^ ||D(/)(m, n) II < 1 — e . 

Now, consider any minimizer w = {wi,W2) with limx^oo'Wi{x) = p+ and limj^^oo 'U^2(— 2;) = P~ ■ 
Then there is an L < 00 so that x > L =^ \wi{x) — p'^\ + \w2{—x) — p^\ < S. Now for x > L + R, 
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Since {wi{x),W2{x)) = <I> * wi{x), * W2{x)j , it follows that ioi x > L + R, \wi{x) — + 

\w2{—x) — < (1 — Iterating this argument leads to the conclusion that for x > L + kR, 
\wi{x) — + \'W2{—x) — p~ \ < {1 — e)^5. A similar argument applies as x tends to — oo. □ 



7 For the functional ^, critical points are minimizers 

We are now ready to prove the main theorem for Q: 

7.1 THEOREM. If {wi,W2) and {vi,V2) are any two critical points of Q in Mp- p+ x M^+ p-, 
then there is an a & R so that 



{vi{x),V2{x)) = {wi{x - a),W2{x - a)) . (7.1) 

Thus, there is exactly one critical point {'Wi,W2) such that wi{0) = W2{0). It is symmetric in the 
sense that wi{x) = W2{—x) for all x. 

Proof: We keep the notation of Section 5. Theorem 15.21 is applicable since by Theorem 16.21 the 
probability densities r]i and i = l,m, have finite moments of every order. Now, if (mx,n\) is the 
displacement convex interpolation between {wi,W2) and {vi,V2), Q{m\,n\) is constant since both 
endpoints are critical points. By the strict convexity up to translation, we see that ()7.ip is true. 

Since {wi,W2) G -A^p-,p+ x -^p+.p-i and both functions are strictly monotonic, there is some b 
such that wi{h) = W2{b). Because of the strict monotonicity of wi; i.e., the strict positivity of its 
derivative, which was proved in Theorem 16.21 this value of b is unique. 

Next, by the symmetries of the functional, since {wi{x),W2{x)) is any minimizer of ^ in A4p- p+ x 
Mp+ p-, then so is {w2i—x),wi{—x)). Hence, by the first part of the Theorem, there is an a G i2 
so that 

{w2{-x),wi{-x)) = {wi{x - a),W2{x - a)) . (7.2) 

Evaluating both sides at x = 0, we see that since u;i(0) = 1^2(0), wi{—a) = W2{—a). By the unique- 
ness of the crossing point established above, a = 0, so that {'W2{—x),wi{—x)) = {wi{x),W2{x)) for 
aU x. □ 



8 Stationary monotone profiles in several dimensions. 

We close the paper by pointing out that our analysis of the two component case can be adapted to 
yield a uniqueness theorem for the one component case in higher dimensions. 

Let Q he a, (d — l)-dimensional cube of size L spanned by the orthogonal vectors ei, . . . ,ed-i 
and Ca^b,n be the set of continuous functions m(x, y) from M x M'^^^ to M such that for all y S R*^"-^ 

lim m{x,y) = a and lim m{x,y) = b , 

X—^ — 00 X^ + OO 

and such that m is L-periodic on R'^^^ in the sense that m{x,y + Lck) = m{x,y) for each k = 
1, . . . ,d — 1 and for each y E M.'^^^. 



CCELM July 2008 



20 



Consider the following d-dimensional free energy on C. 



— a,a,S7 



!F{m) = / F{m{x,y))<lxdy 



(8.1) 




RxQ jRxQ 



(m(xi,yi) - m{x2,y2)fj{xi - X2,yi - y2)d2;idx2dyidy2 , 



J(^x,y) = U{^J x"^ + with U monotone decreasing, finite range smooth function on [0, +00) 

and F an even double well potential with minima in —a and a and F(iba) = 0. (These specific 
conditions on F enable us to be brief, and can easily be relaxed.) 

Obviously, if fh{x) is a minimizer for the corresponding one dimensional problem, then 



is a critical point of J- on C_a,a,n, and is an obvious candidate to be the unique minimizer. We shall 
show here that not only is it the minimizer - this fact has been proved by Alberti [1] - but that, 
up to translation in x, fh{x, y) is the only solution of the Euler-Lagrange equation for minimization 
of JF that is monotone in x for all y. A related question as to whether all monotone solutions of 
the Euler-Lagrange equation have this special form has been extensively investigated for the local 
variant of the free energy (Allen-Cahn or van der Waals) with J \\/m{x, y) p in place of the non-local 
interaction integral above. It turns out that the non-local case may be easily treated by regarding 
the one dimensional profiles x 1— > m{x, y) for different y as profiles for different components, and 
applying our previous results. 

Define Ma,b,n to be the subset of Ca^b,n for which m{x,y) is monotone in x for each y € M"'^^. 
As before, there is a rearrangement inequality that allows one to reduce the minimization problem 
over C-a^a,n to minimization over ^A-a,a,n■ Given m € C_a,a,r2 we define m* € ^A-a,a,n as follows 
by separately rearranging m{-,y) for each y £ M"^"^, using the one dimensional rearrangement 
procedure. By the rearrangement results cited above, J-{m*) < F{m). In [1], Alberti proceeds 
with a careful study of the cases of equality here. Instead, we henceforth restrict our attention 
to m € M-a,a,Q.^ and shall show that up to translation in x, there is just one solution of the 
Euler-Lagrange equation in this set. 



where J = /^xr "^(^j y)dxAy. We now observe that the second term on the right can be written in 
terms of the I functional that has been studied in Section [5j Indeed, this term can be written as 



m(x, y) := m{x) 







This identity relates the multidimensional problem to the two species problem: here m( • , yi) plays 
the role of the profile for one species, and —m{-,y2) plays the role of the profile for the other 
species. 
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Now notice that for a + b = (or c + d = 0), the statement of Lemma 15.11 simpHfies in a 
significant way: The first moments drop out as in Corollary I5.3| and we have (using the notation 
from the lemma) 

-J(mi, -1712) =Aa^j / W{xi - X2)pi{xi,yi)p2{x2, y2))dxidx2 + Qa^a . 

(First moments could be dealt with as before, but we avoid doing so in order to focus on how one 
may regard the multidimensional problem as a multi-component problem, which is the main point 
of this section.) 

Given two profiles rriQ and mi in M.-a,a,n, let m\ be the interpolation defined by interpolating 
between mo( • , y) and mi( • , y) separately in each y. Let x i-^ T{x, y) be the corresponding optimal 
transportation plan, and let S{x,y) = T{x,y) — x. Let m\{x,y) be the induced interpolation 
between mQ{x,y) and mi{x,y). Then 

- 'I{mx{ ■ , yi), -mx{ ■ , 2/2)) = 

4o^ / / W[xi -X2 + X{S{xi,yi) - S{x2, y2)]pi{xi,yi)p2{x2, y2)dxidx2 + 6a^a . 

JR JR 

Since W is strictly convex near the origin if J is strictly positive near the origin, it follows that 
if y2 and yi sufficiently close to one another, pi{xi,yi)dxi and /52(a^2; y2)dx2 both assign positive 
mass to some small interval around some xq. Therefore, for such yi and y2, we see that A 1— > 
—Z{mx{ ■ ,yi), —TTixi ■ 72/2)) is strictly convex, and for any yi and y2 it is convex. Clearly, the set 
of points (^1,2/2) for which we have strict convexity is a set of positive measure (containing the 
diagonal) with respect to dyidy2, and so 



A^-/ I{mxi- ,yi),-mx{- ,y2))dyidy2 
Jnxn 

is strictly convex, apart from translation in x. This strict convexity proves that, up to translation 
in X, there is just one critical point of J~ in ■AA~a,a 

^Q. Since clearly m{x,y) is a critical point, we 

have the following: 

8.1 THEOREM. Assume that J is bounded below by a strictly positive number on some neighbor- 
hood of the origin. Let m{x, y) be any solution of the Euler- Lagrange equation for the minimization 
of T that belongs to M.-a,a,n- Then for some xq € M, m{x,y) = fh{x — xq) for al x and y, where 
fh the antisymmetric minimizer for one dimension. 
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